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Abstract 


A  detailed  account  o£  the  research  accomplished  as  per  the 
-tatement  of  or^  of  this  conLract  is  i_,iven  herein.  Of  particular  inter¬ 
est  to  the  needs  of  Project  VEL4-UNIFf‘RiV.  are  the  folloA'ing;  a  study 
of  the  radiation  patterns  of  surface  waves  from  underground  nuclear 
explosions  and  siTall  magnitude  earthquakes,  wherein  it  was  determined 
that  the  radiation  patterns  frorr!  explosions  in  tuff  and  alluvium,  hut  not 
in  granite,  could  be  interpreted  as  arising  from  an  explosive  force 
acting  as  a  step  function  in  time;  a  study  of  the  relative  excitation  of 
surface  waves  by  earthquakes  and  underground  explosions,  which 
yielded  the  result  that  most  but  not  all  of  the  earthquakes  studied 
generated  surface  waves  5  to  10  times  greater  than  the  rraximurn  ob¬ 
served  for  explosions  of  the  same  Richter  miagnitude;  an  almost  com¬ 
plete  summary  of  Rayleigh  wave  phase  and  group  velocities  for  various 
regions  of  tht  world;  a  programi  of  epicenter  relocation  carried  out  so 
far  for  the  South  Pacific  Ocean,  the  Arctic  north  of  bO^^N.  the  Gulf  of 
-.den,  Red  Sea,  and  Arabian  Sea,  the  African  Rift  /  one,  and  the  Central 
Indian  Ccean,  which  shoved  that  in  these  places  at  least  large  earth¬ 
quakes  are  confined  to  much  smaller  regions  than  previously  thought; 
successful  utilization  of  displacement  transducers  with  long  period 
seismometers,  thereby  improving  the  long  period  response  of  these 
instruments  by  removing  the  6  db/octave  fall-off  v/ith  increasing  period, 
which  is  inherent  in  velocity  transducers. 
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At  present,  the  current  Statement  of  ^Jork  for  this  contract 
reads  as  follows: 

The  Contractor  shall  sunply  the  necessary  personnel,  facilities, 
servic-^s  and  materials  to  accomplish  the  follovdng; 

Item  1  -  Extension  of  e>d.stin^  seismograph  pier  to  accommodate 
the  more  sensitive,  longer  period  seismological  instrumentation  required 
to  carry  out  the  research  herein. 

Item  2  -  Investigation  of  the  applicability  to  intermediate 
(2  to  10  sec)  and  long  (10  sec  and  greater)  period  seismic  leaves  of  phase 
compensation  techniques  for  the  purpose  of  detemir-ing  an  effective 
source  polarity  as  a  function  of  azimuth.  In  this  connection  the  com¬ 
pilation  of  world-wide  charts  or  tables  of  phase  and  group  velocity  for 
surface  waves  of  various  periods  will  be  initiated. 

Item  3  -  Development  of  high-speed  computational  techniques 
for  the  solution  of  the  normal-mode  propagation  problems  and  the  Fourier 
analysis  and  synthesis  problems  that  arise  in  connection  with  Item  2, 

Item  U  -  Development  of  seismic  instruirentation  with  increased 
sernitivity  in  the  long  period  ranee  ard  with  improved  discrimination 
against  noise  in  the  dominant  microseism  band. 

Item  5  -  Conduct  studies  and  develop  techniques  in  the  use  of 
specialized  instrumentation  for  seismic  data  recording  and  analysis. 
Included  would  be  investigations  of  tape  recording  systems,  filtering, 
digital  analysis  and  analogue  computers  for  the  purposes  of  improving 
detection  techniques  and  methods  of  analysis  of  seismic  records. 
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Item  6  -  Conduct  tests  to  increase  the  knovjledge  of  crustal 
structure  by  land  refraction  iiethods.  Field  studies,  using  a  portable 
seismic  recording  system,  xd.ll  be  conducted  at  selected  sites. 

Item  7  -  At  a  site  to  be  designated  by  the  Contracting  Officer, 
establish  a  long  period  seismograph  station  to  form,  with  the  exi^  ing 
Palisades  and  Uaynesburg  Stations,  a  tripartite  array  for  basic  studies 
of  bot^  and  sxirface  waves.  An  additional  station  within  the  large  array 
m?ty  be  established  to  form  a  smaller  inner  array  if  these  stvdies  indicate 
it  to  be  advantageous. 

Item  0  -  Develop  techniques  for  analyzing  microseismic  data  and 
conduct  studies  of  the  microseismic  spectra* 

Note  1  -  In  the  performance  of  tte  above  research,  the  contractor 
vTill  continue  to  operate  the  widely  distributed  network  of  long  and  inter¬ 
mediate  period  seismic  stations  that  were  initially  put  into  operation  in 
connection  with  the  International  Geophysical  Year,  with  such  changes  in 
numbers,  location,  or  instrumentation  as  may  be  deemed  desirable  for  the 
performance  of  the  contract. 

The  work  accomplished  uxider  each  item  x^ill  be  described  in  detail 
in  the  bulk  of  this  report. 

Item  1  -  Extension  of  existing  seisirragraph  pier  to  accommodate  the  more 
sensitive,  longer  period  seismological  instrumentation  required  to  carry 
out  the  research  herein. 

Bids  xrere  obtained  for  extending  the  seismograph  pier,  and 
also  for  building  a  amall  additional  pier  at  another  location  on  the 
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Lanont  grounds.  It  Tfss  doternined  th'’,t  the  second  alternative  was  far 
less  costly,  and  accordingly  the  Siiall  additional  pier  was  constructed, 
•and  is  in  use  for  long  period  seismographs. 

Item  5  -  Investigation  of  the  applicability  to  intermediate  (2  to  10  sec) 
and  long  (10  sec  and  greater)  period  seismic  waves  of  phase  compensation 
techniques  for  the  purpose  of  determining  an  effective  source  polarity  as 
a  function  of  azimuth.  In  this  connection  the  compilation  of  vrorld-wide 
charts  or  tables  of  phase  and  group  velocity  for  surface  xmves  of  various 
periods  will  be  initiated, 

A,  Radiation  Rxtterns 

A  paper,  entitled  "Radiaticn  Patterns  of  Surface  ’Javes  from 
Undergound  I'luclear  Explosions  and  Small  Magnitude  Earthquakes"  by 
James  N,  Brurje  and  Paul  W,  Pomeroy,  v.as  submitted  to  the  Journal  of 
Geophysical  Research, 

In  this  paper,  the  radiation  patterns  of  surface  waves  were 
used  to  study  the  source  mechanisms  of  undergrouri  nuclear  explosions, 
associated  collapses  and  small  magnitude  earthquakes.  This  study  is 
the  first  successful  attempt  to  use  the  technique  of  Fourier  analysis 
and  phase  equalization  as  a  function  of  azimuth  to  study  the  source 
mechanisms  of  underground  nuclear  explosions.  The  results  are  xianortant 
to  the  problem  of  distinguishing  urderground  nuclear  explosions  from 
earthquakes  on  the  basis  of  seasmiu  wave  generation  as  well  as  to  the 
problem  of  determining  earthquake  source  mechanisms.  The  successful 
application  of  the  technioues  used  in  this  paper  for  the  study  of 
explosion  source  mechanism,  uliich  in  many  vrys  is  simpler  to  study 
than  for  earthquakes,  indicates  that  with  iiaproved  epicentral  information, 
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better  instruir^ntation  and  accurate  regional  phase  velocity  data, 
these  techniques  uill  ultimately  yield  extremely  valuable  information 
about  the  source  ireclianisms  for  earthquakes  ard  explosions  and  about 
the  tectonic  strains  within  the  earth. 

In  recent  years,  several  authors  have  used  surface  waves  to 
study  source  mechanism.  However,  because  of  a  lack  of  a  sufficient 
number  of  adequately  equipped  seismopTrajii  stations,  the  method  has  not 
been  fully  utilized.  The  temporary  Ii?SM  (Long  Range  Seismic  Measure¬ 
ments)  stations,  established  to  monitor  the  xmderground  nuclear  ex¬ 
plosions  at  the  Nevada  Test  Site,  operated  long  period  instruments 
sharply  peaked  around  25  seconds,  ''ud  often  clearly  recorded  the 
fundamental  mode  of  Love  and  Rayleigh  waves  in  the  period  range  from 
10  to  30  seconds.  Thus  they  provide  an  excellent  source  of  data  for 
the  study  of  the  source  irechcnisms  of  underground  nuclear  explosions, 
since  they  surround  the  tpst  site  and  f-e  horizontal  instruments  are 
oriented  in  radial  and  transverse  directions  to  the  test  site. 

Seismograms  from  these  stations  for  the  larger  underground 
nuclear  explosions  and  certain  earthquakes  vrere  studied.  Selected 
records  were  digitized  and  the  digital  records  were  then  Fourier 
analysed  to  determine  amplitude  and  phase  spectra.  The  phase  spectra 
were  then  equalized  to  tae  source.  The  results  of  this  analysis  were 
studied  as  a  function  of  azimuth  for  each  event.  The  results  were  then 
interpreted  in  terms  of  possible  source  mechanisms. 

It  was  concluded  that: 

1,  The  radiation  patterns  of  Rayleigh  waves  generated  by 
underground  nuclear  explosions  in  tuff  and  alluvium  suggest  that  the 
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source  mechanism  is  an  explosive  force  acting  as  a  step  function  in  time. 
The  time  function  agrees  with  that  obtained  ty  other  authors  from 
doubly  integrated  accelerometer  recoids,  Hovrever,  the  existence  of 
v/aves  and  small  asyjnmetries  in  Rayleigh  vmve  amplitudes  indicate 
that  more  complicated  forces  also  act  at  the  source  for  these  events, 

2,  The  collapses  which  follow  many  of  the  larger  explosions 
in  tuff  and  alluvium  generate  Rayleigh  waves  with  reversed  polarity 
relative  to  the  explosions.  The  source  mechanism  could  be  an  implosive 
force  acting  as  a  step  function  in  time  (or  a  dotjmrard  force  acting  as 
a  step  function  in  time).  Love  wave  generation  by  the  collapse  events 
is  ejctrem^ly  small  compared  to  that  by  the  explosions.  This  fact 
indj.cates  that  the  generation  of  long  period  Love  waves  is  a  property  of 
the  source  for  the  explosions,  and  is  not  a  result  of  mode  conversion 
along  the  propagation  path, 

3.  The  explosion  in  granite  (Hardhat)  produced  a  surface 
wave  radiation  pattern  with  nearly  perfect  double-couple  sjammetry  in 
both  phase  and  amplitude.  Such  a  pattern  could  not  be  generated  by  a 
simple  explosive  force  and  must  have  been  caused  ty  cracking,  forced 
motion  along  pre-existing  faults,  or  by  tectonic  strain  release  due 
either  to  release  of  strain  in  the  volume  of  rock  shattered  by  the  ex¬ 
plosion  or  triggering  of  tectonic  miOtion  along  nearby  faults. 

ij.  The  amplitude  and  initial  phase  spectrums  of  Rayleigh 
waves  from  earthquakes  a.re  usually  more  complicated  than  the  spectrums 
of  Rayleigh  waves  from  the  explosions  in  tuff  and  alluvium.  The 
radiation  pattern  of  surface  vraves  for  the  .irisona-l?tah  border  earth¬ 
quake  of  1^  February  1962  showed  an  initial  phase  of  about  zero  for 
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Rayleigh  waves  at  all  asL'miths  for  which  it  was  possible  to  make 
reliabla  r^asureiaeni^.  This  is  opposite  in  phase  to  the  phase  ob¬ 
served  for  the  explosions  in  tuff  and  alluvium.  It  is  expected  that 
the  radiation  patterns  for  Rayleigh  waves  generated  hy  eartiiouakos  will 
be  quite  different  from  the  radiation  patterns  for  explosions  in  tuff 
and  allmdum  and  the  results  of  this  iitud;^'  indicate  that  in  the  future 
it  will  be  possible  to  use  the  inethod  described  here  to  stidy  the  source 
mechanisms  of  numerous  eai'thquakes . 

B.  Si^ctral  Analysis  of  Body  Vlaves  and  Synthesis  of  Time 
Function 

A  method  to  compute  the  excitation  time  function  at  the  origin 
has  been  developed.  This  nsthod  consists  of  the  follotong  processes? 

1.  The  instrumental  characteristics  of  an  electromagnetic - 
seismofuraph  have  been  computed  and  the  effect  of  the  instrument  may  be 
removed  froru  a  spectrum  of  a  recorded  body  wave.  A  true  ground  motion 
is  thus  obtained, 

2,  It  is  knovm  that  the  spectrum  of  the  ground  mntion  of  a 
bocfy  wave  consists  of  tv70  parts,  a  spatial  spectrum  S(p)  and  a  tim^ 
spectrum  F{p),  A  mathematical  expression  for  the  spatial  spectrum  may 
be  derived  for  a  spherical  origin. 

By  removing  the  spatial  spectrum  which  is  computed  uixler  an 
assumption  of  a  double-couple  type  stress  distiibution  from  the  spectrum, 
ihe  time  spectrum  i(p)  can  bo  obtained  from  a  spectrum  of  an  observed 
ground  moveiiient.  Finally,  the  tima  function  is  synthesised  by  maans  of 
nuirierical  integration. 

Synthesised  tiiae  functions  have  been  c Diluted  for  12  earth- 


quakes  (i>  shallow  earthquakes,  8  dcep-focus  earthquakes).  It  is  fou 
that  the  tine  functions  of  deep-focus  earthquakes  begin  with  a  sharp 
conmencerr^nt  and  die  out  very  gradually,  or  look  like  step  functions 
eicpected.  The  tine  functions  of  shallo’.?  earthquakes  are  generally  m 
complicated!  most  of  them  shov?  a  doublet  type  of  variation  -and  some 
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them  are  oven  more  oscillatory. 


This  difference  of  the  time  fwictions  between  shallov?  and  deep- 
focus  earthquakes  implies  that  the  effect  of  the  free  surface  has  an 
important  role  for  the  shallow;  earthquakes.  To  study  this  problem,  the 
reflection  of  the  waves  caused  by  a  free  surface  near  the  orig,.n  is 
talffln  into  account  for  the  computation  of  a  t*.me  function  of  a  shallow- 
focus  earthquake.  It  is  found  that  some  of  the  time  functions  of  shallow 
earthquakes  become  very  similar  to  those  of  deep-focus  earthquakes  vjhen 
the  reflection  caused  by  a  frC'.:.  surface  is  removed* 


Fl'osumabl^'  we  may  conclude  that  there  is  no  essential  difference 


in  the  mechanism  of  the  strain  release  for  shallow  and  for  deop-focus 


earthquakes  except  for  some  major  shallow  earthquakes  where  it  was  im¬ 
possible  to  obtain  a  simple  fin-e  function  even  if  a  free  surface  re- 


’uction  lias  been  dene. 


Also,  an  attempt  to  compute  a  time  function  of  an  eOTlosion 
has  been  made  {Hardhat,  February'-  15,  1962).  The  s^withcsiced  time  function 
is  also  not  simple  in  this  case ,  It  is  believed  that  the  frequency  range 
in  this  analysis  is  higher  than  for  the  previous  analyses  and,  accordingly, 
minor  structure  of  the  crust  has  an  immurtent  effect  on  the  free  surface 
reduction.  It  is  difficult  to  make  an  acdorate  discussion  of  the  excit¬ 
ation  function  of  an  explosion  unless  we  know  an  ex-act  free  surface 


reduction. 
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u,  Qsneral  Smanury  of  Phase  Velocity  flarninp  rYojr'st 

Rayleigh  vjave  phase  velocities  are  generally  determined  1^ 
the  f  olloidng  iKithods  : 

1»  Tripartite-net  method  by  correlating  the  same  phase  in  a 
three-station  net* 

2,  Two-station  correlation  method  by  selecting  at  least 
two  stations  which  lie  approximately  in  a  common  great  circle  v/ith  the 
epicenter.  The  wave  front  is  assuiiicd  to  be  perpeniicular  to  the  great 
cirle  path* 

3*  One-station  m^ethod  by  assuming  the  initial  phase  to  be 
indepemient  of  the  frequency  at  the  origin. 

Phase  velocities  thus  determined  by  marXr'  invest! -^ators  in 
various  parts  of  the  world  provide  the  principal  data  for  the  present 
report.  Haps  have  been  prepared  in  the  forms  of  regional  velocity  maps 
and  velocity-contour  raps, 

a,  Eurasia-Africa 

The  phase  velocity  values,  which  we  have  obtained  for  this 
region,  are  determined  by  means  of  the  two-station  correlation  and  the 
one-station  methods.  They  have  long  paths  across  several  geological  or 
physio  graphical  provinces,  Hov/ever,  the  phase  velocities  through  Tibet 
and  Himalaya  are  lower  than  those  for  the  oaths  across  the  nerthern  part 
of  USSR,  The  phase  velocities  for  Africa  correscorxi  to  standard  conti¬ 
nental  structure,  i.e,  about  35  km  in  the  crustal  tMckness.  These 
average  phase  velocities  give  only  a  gross  picture  of  the  average 
structure  for  the  wave  propagation  paths. 
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Phase  velocities  have  heen  laeasorGd  by  th^  tripartite -net 
iTiothod  in  the  Alps  aM  Mediterranean  by  Plnopoff  and  Press*  Their  data, 
wbJ-Ch  were  determined  chiefly  from  Atlantic  earthquakes,  are  used  in 
our  regional  velocity  map  as  these  waves  show  less  interference* 

b»  Pacific  Ocean 

Kuo,  BruCie  arA  Major  used  the  one-station  and  tuo-station 
correlation  methods  to  determine  phase  velocities  for  fourteen  paths 
across  the  Pacific  centered  at  Suva.  Pcr»roy  determLned  the  ph^se 
velocity  between  Bikini  and  Eniwetok.  Ra;i/leigh  nave  phase  velocities 
for  the  paths  through  the  northern  Pacific  in  the  deep  ocean  basins  are 
definitely  higher  than  for  those  patlis  across  the  southeastern  Pacific, 
vfnere  the  Easter  Island  Rise  is  situated.  Phase  velocities  are  much 
lower  for  those  paths  across  the  Ifelanesian  arA  Tasman-ffew  ZealaTid 
tectonic  complex, 

c.  Hew  Zealand 

Thomson  and  Evison  published  their  phase  velocity  data  for 
North  Island  and  the  northern  part  of  South  Island  in  the  neried  range 
from  20  to  30  seconds.  Their  preliminary  results  seem  to  indicate  that 
phase  velocities  are  quite  uniform  in  this  cart  of  Hew  Zealand, 

d.  Japan 

Aki’s  phase  velocity  deterrdnation  by  the  tripartite -net 
method  in  Japan  for  the  period  of  29-uecond  rtaylcich  waves  gives  a  good 
correlation  betvreen  topography  and  phase  velocity.  Lew  phase  velocity 
coincides  ndth  the  htigher  elevation  in  Japan. 

e.  South  Asarica 

This  study  is  still  in  pro-ress  at  lament.  Phase  velocities 
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in  South  Aiaerica  are  chiefly  determined  by  the  txTO-station  correlation 
method.  The  preliirinaiy  phase  velocity  data  is  consistent  with  respect 
to  the  geological  provinces.  The  paths  containing  a  large  percentage 
in  the  Andes  have  lower  phase  velocities  than  those  containing  a  small 
percentage  in  the  Andes , 

f .  Eastern  Canada 

Brune  aid  Dorman  found  phase  velocities  are  relatively  high 
in  the  Canadian  shield  area, 

S*  U,  A, 

The  irnst  complete  data  available,  at  present,  is  in  the  U.  S.  A. 
Ewing  and  Press  in  1959  determined  phase  velocities  across  the  United 
States  by  using  a  distant  earthquake  of  Samoa  in  the  Pacific.  Th.eir 
method  was  chiefly  tl»  tripartite -net  method.  Their  work  has  been  ex¬ 
tended  here  to  include  the  most  recent  data  of  the  lESM  stations.  In 
the  present  extension  an  Atlantic  earthquake  was  used  so  that  a  reversed 
Rayleigh  wave  propagation  pattern  through  the  U,  S,  could  be  established. 
To  determine  the  phase  velocities  in  the  Appalachian  l-fountain  area,  the 
Mexican  earthquake  of  September  1^2  was  used,  since  th^re  were  six 
IRSM  stations  in  operation  at  that  tin®  and  these  stations  are  approxi¬ 
mately  in  a  common  great  circle  with  the  epicenter.  In  th;  Pennsylvania 
and  New  York  area,  phase  velocities  were  determimd  by  Oliver,  Idvach 
and  Dorman. 

Cur  most  rcct^nt  results  of  phase  velocities  determined  by 
both  the  tripartite -net  ard  the  two-station  correlation  methods,  arc 
in  good  agreement  with  these  previous  results  published  by  Swing  ani 
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Press  Oliver^  Kovach  and  Dorman.  Tne  phase  velocities  are  highest 
at  the  West  Coast  and  lowest  in  the  Basin  and  Range  Province,  Colorado 
Plateau,  and  Rocky  Mowitains,  The  phase  velocities  become  higher, 
in  comparison  with  the  high  land  toward  the  cast  in  the  Interior  Plain, 
and  again  decrease  along  the  Appalachian  i^imtains. 

With  these  data,  velocity-contour  maps  of  the  ITiited  States 
were  constructed  at  5-second  intervals  for  20,  25,  30,  and  35-second 
Raylei^  vjaves. 

Detailed  studies  of  the  crust  and  upper  mantle  structure  of 
the  United  States  were  made  along  two  profiles  across  the  continent  on 
the  basis  of  the  phase  velocity  data,  seismic  refraction  data  and 
gravity.  These  te-Jo  profiles  are  (1)  from  Southern  California  tijough 
Arizona,  New  Mexico  and  Tennessee  to  the  Southern  Appalachians,  and 
(2)  at  the  northern  port  of  the  United  States  extending  from  San 
Francisco  through  the  Basin  and  Range  Rrovince,  Rocky  Mountains,  and 
the  Interior  Plains  to  New  York. 

ITie  phase  velocity  data  support  the  existence  of  the  inter- 
irjediate  layer  in  the  crust,  which  is  in  general  identified  by  second 
arrivals  in  the  seismic  refraction  method.  Nevertheless,  for  the 
Pennsylvania  and  Ifew  York  area  and  for  the  west  coast  of  California, 
the  ass'omed  one-layer  structure  also  fits  the  observed  data.  This 
demonstrates  that  the  structure  interpretation  is  never  unique. 
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IteiTi  3  -  Development  of  high-speed  coirpsutational  techniques  for  the 
solutj.on  of  the  normal -mode  propagation  problems  and  the  Fourier 
analysis  and  synthesis  problems  that  arise  in  connection  with  Item  2. 

A.  Travel  Times,  Body  Waves  and  Dispersion  in  the  Earth 
A  simple  relation  between  body  waves  and  higher  mode  dispersion 
has  been  derived,  based  on  a  constructive  phase  interference  requirement 
for  body  waves  such  as  S  and  PS.  For  a  normal  .mode,  the  travel  time,  tt, 
for  a  given  distance,  ^  and  phase  velocity,  l/(dt/dZ\),  are  related 
to  the  period,  T,  by  the  equation: 

tt  -  A  (dt/d  A)  “  n  /2  7/  )T 

where  n  is  an  integer  equal  to  the  mode  number  minus  one,  and  *5  is  a 

phase  shift  arising  from  reflection  and  propagation  of  the  wave.  The 

group  velocity  is  given  by  /(tt  -  (g*^)  »  const  .  For  a  homogeneous 

sphere  4  is  the  phase  shift  for  reflection  from  a  spherical  surface. 

^  (and  T)  may  be  experimentally  determined  from  the  difference  in  phase 

of  Fourier  transforms  of  successive  surface  reflections  of  body  waves 

corresponding  to  rays  of  a  given  phase  velocity  (e.g,  ScS  at  one  station 

and  ScSScS  at  another  station  at  txdce  the  distance).  The  method  may  be 

used  to  measure  torsional  dispersion  from  S„  and  S  cS  waves  ani  over 

h  H  H  * 

a  limited  phase  velocily  range,  spheroidal  dispersion  from  IS  waves. 

The  discovery  explains  very  simply  the  relationship  between 
body  waves  and  higher  normal  nodes  and  thus  adds  greatly  to  the  under¬ 
standing  of  wave  propagation  in  the  earth* 

Experimental  values  of  dispersion  are  given  for  several 
values  of  phase  ^^elocity  (8  to  3^  km/sec)  am  modes  from  3  to  more  than 
10  and  compared  with  theoretical  values. 
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3 ,  Calculation  of  Free  Periods  of  Oscillation 

Two  programs  for  the  IBM  7090  have  been  completed  for  the 
study  of  the  free  oscillation  of  the  earth.  One  oronram  computes  the 
periods  of  the  free  vibrations  for  earth  models  consisting  of  a  solid 
mantle  and  a  liquid  core.  The  other  program  has  provision  for  treating 
earth  models  with  a  solid  inner  core.  In  additionj  both  programs 
compute  phase  and  group  velocities  and  the  amplitudes  of  radial  and 
trarsverse  particle  motions  and  the  gravitational  perturbation  as  a 
function  of  depth. 

The  free  periods  and  resulting  phase  velocities  have  been 
calculated  for  all  orders  vhlth  free  periods  greater  than  2C0  seconds 
for  the  fuTidamental  and  first  ard  secorai  higher  modes  for  four  earth 
models.  The  four  models  used  vrere  obtained  by  combining  compressional 
and  shear  velocities  according  to  Jeffreys  or  according  to  Gutenberg 
with  either  Eullen’s  density  model  A  or  B.  The  Bullen  B  model  used 
differed  from  that  originally  proposed  ty  Bullen  in  that  the  inner  core 
was  assumed  to  to  liquid,' 

V/ith  these  data  it  is  possible  to  obtain  theoretical  phase 
velocity  curves  for  long  period  surface  and  mantle  xfaves  above  2CX) 
secoMs.  Uith  increasing  period  these  waves  become  less  and  less  de¬ 
pendent  on  crustal  variations  and  therefore  it  is  much  easier  to  deter¬ 
mine  the  source  mechanism,  Lhfortunately,  in  this  case  theory  is  ahead 
of  observational  data  and  must  wait  for  more  sensitive  long  period 
seismographs. 

The  effect  on  the  free  periods  of  introducing  rigidity  into 
the  inner  core  has  also  been  studied.  The  core  modes  are  very  sensitive 


to  the  rigidity  of  the  inner  core,  decreasing  rapidly  in  period  with 
increasing  shear  velocity  of  the  inner  core.  However,  there  is  very 
little  effect  on  the  norriSl  types  of  spheroidal  vibrations,  except 
v/hen  the  period  is  close  to  that  of  a  core  mode.  Under  these  circum¬ 
stances  an  interaction  occurs  between  the  two  modes  in  which  the 
vibrations  take  on  characteristics  intermediate  between  the  two  modes, 
and  an  appreciable  shift  in  period  occurs. 

During  the  course  of  this  investigation,  several  new  modes 
of  vibration  have  been  discovered.  One  of  these,  a  long  period  vibration 
of  order  number  one,  is  very  similar  to  the  vibration  propoeed  by 
Slichter  to  explain  the  period  of  86  minutes  obtained  in  the  spectrum 
of  free  vibrations  from  the  gravimeter  record  of  the  great  Chilean 
earthquake  of  I960, 

C,  Rayleigh  Wave  Phase  Velocities  at  Long  Periods 

Spectral  peaks  corresponding  to  the  free  periods  of  oscillation 
of  the  earth  have  been  found  in  the  spectra  of  eight  seismograms  obtained 
vertical  seismographs  locate'^  at  various  stations  around  the  vjorld. 
The  stations  are  all  part  of  the  Lament  world-wide  netv/ork  of  long- 
period  seismographs  initially  placed  in  operation  during  the  International 
Geophysical  Year,  The  stations  in  question  are  located  at  Agra,  Indiaj 
Hallett,  Antarcticaj  Hong  Kong;  Lwiro,  Republic  of  the  Congo;  Mb, 

Tsukuba,  Japan;  Resolute,  Korthern  Canadian  Islands;  Suva,  Fiji;  and 
Uppsala,  Sweden.  All  of  the  data  were  obtained  following  the  great 
Chilean  earthquake  of  22  May  I960,  The  vertical  seismographs  are  all 
equipped  with  l5-second  pendulums  aixi  75-3econd  galvanomsters , 


The  records  used  started  from  approxijnately  four  to  approxi- 
iTiately  fifteen  hours  after  th^  earthquake  occurred.  In  most  cases  the 


seismometers  were  driven  off  scale  and  were  reset  the  following  day 
when  the  records  were  changed.  This  is  the  primary  reason  for  the 


different  beginning  times  of  the  records.  As  much  as  possible  of  one 
complete  day *3  record  was  used  for  the  analysis, 

Tlie  seismograms  were  digitised  at  Lanont,  All  of  the  records 


were  digitized  with  an  interval  of  approximately  six  secords  befe’een 
data  points.  The  exact  value  of  the  digitisation  interval  was  cal¬ 
culated  to  three  decimal  places  by  dividing  the  time  length  of  the 
record  by  the  number  of  data  points.  In  tbds  way  allowance  could  be 


made  for  the  stretching  or  shrinking  of  the  record  paper  during  processing. 

Prior  to  digitising,  the  records  were  smoothed  by  eye  to  re¬ 
move  extremely  short  period  waves.  After  digitisation,  the  data  were 
smoothed  by  averaging  over  19  and  21  points,  and  then  every  third  point 
was  selected  in  order  to  reduce  the  amount  of  the  data  while  avoiding 
aliasing.  The  data  were  then  Fourier  analyzed  and  the  periods  of  the 
spectral  peaks  obtained  by  interpolation. 

The  various  spectra  show  the  features  to  be  expected.  Since 
these  are  all  the  spectra  of  vertical  seismograms,  cnly  the  soheroidal 
types  of  oscillation  are  observed.  The  general  appearar»ce  of  each  of  the 
spectra  is  approximately  the  same.  The  rapid  fall-off  of  the  sensitivity 
of  the  seismograph  at  long  periods  limits  the  longest  observed  free 


oscillations  to  those  trlth  periods  of  about  600  to  7C0  seconds,  that  is, 
to  approximately  orders  8  and  9  of  the  fundamental  mode.  The  shortest 
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period  oscillations  observed  are  of  the  order  of  about  200  seconds. 

As  might  be  expected,  the  short  period  oscillations  decay  more  rapidly 
\^ith  tirse  than  the  long  period  oscillations.  Therefore,  the  spectra  of 
the  seismograms  starting  many  hours  after  the  earthquake  do  not  have 
peaks  with  as  short  a  period  as  the  spectra  of  seismocrams  starting  a 
few  hours  after  the  earthquake. 

The  highest  order  observed  for  the  first  higher  n^de  is 
possibly  order  9*  This  is  In  accord  with  theoretical  calculations  which 
indicate  that  the  amplitudes  of  orders  9  through  l5  will  be  small  at  the 
surface  of  the  earth. 

The  data  discussed  so  far  have  been  combined  with  all  other 
free  oscillation  data  available  for  the  spheroidal  oscillations  to 
obtain  average  values  for  the  periods  of  the  free  oscillations.  There 
is  a  slight  problem  in  choosing  data  about  spheroidal  oscillations. 

For  orders  10  and  below  of  the  fundamental  mode  there  can  be  no  confusion 
in  distinguishing  spheroidal  from  torsional  oscillations  with  the  ex¬ 
ception  of  and  so  that  all  of  the  data  available  were  used  for 
these  orders.  For  orders  higher  than  10  it  becomes  difficult  to 
positively  identify  a  spectral  peak  as  belonging  to  a  spheroidal  or  to 
a  torsional  oscillation.  Accordingly,  for  orders  higher  than  10,  only 
data  obtained  by  vertical  seismographs  or  by  gravimeters,  both  of  which 
record  only  spheroidal  oscillations,  were  used. 

The  data  selected  in  accordance  with  the  above  criteria  In 
addition  to  that  pa'eviously  discussed  x^ere  obtained  ty  workers  at  Los 
Angeles j  Pasadenaj  Palisades j  Chester,  New  Jersey:  Tiefenort,  Germany j 
Trieste}  Paris}  and  lycto.  The  data  of  tlte  first  four  were  used  for 


1? 


all  orders,  of  the  last  four  only  for  orders  10  and  below. 

Since  the  earth  is  neither  homogeneous  nor  a  sphere,  arid  its 
elastic  parameters  and  density  vary  azimithally  as  well  as  radially, 
it  is  not  to  be  expected  that  the  periods  of  free  oscillation  will  be 
the  same  at  all  stations.  This  may  be  seen  intoitivel;y  by  considering 
the  free  oscillations  as  standing  waves  which  arise  from  the  constnictive 
interference  of  surface  waves.  Thus,  if  the  great  circle  paths  through 
the  several  stations  differ,  it  can  be  expected  that  different  periods 
of  surface  waves  ^jill  constructively  interfere  among  the  different  paths. 
Therefore,  a  systematic  variation  between  stations  should  be  expected, 
particularly  for  the  shorter  period  oscillations  xvhich  are  confined  to 
the  crust  and  upper  mantle.  A  comparison  of  the  observed  periods  for 
the  different  stations  with  the  average  period  does  show  this  systematic 
difference  at  the  shorter  periods,  that  is,  200  to  350  seconds. 

The  phase  velocities  of  Rayleigh  waves  can  be  calculated  from 
the  periods  of  the  spheroidal  oscillations.  The  phase  velocities  of 
course  can  be  calculated  only  at  the  periods  of  the  free  oscillations. 
The  actual  curve,  continuous  from  200  to  3200  seconds,  is  obtained  by 
interpolating  graphically.  Over  raich  of  the  range  the  standard  error 
of  the  mean  is  less  than  one  part  in  a  thousand. 

For  purposes  of  comparison,  phase  velocity  curves  obtained 
from  the  calculated  free  periods  of  oscillation  for  two  earth  models 
have  been  compared  with  the  observed  data.  The  models  both  have 
velocities  according  to  Gutenberg  but  one  has  ths  densities  of  BuUen’s 
model  A,  the  other  of  Bullen's  model  B.  For  periods  greater  thar.  ICOO 
seconds,  the  phase  velocities  calculated  from  the  observed  free  periods 


agr  rery  well  with  those  calculated  for  tte  Outenberg-Bullen  B  model. 
Per  pt.icds  less  than  ICOO  seconds,  the  observed  phase  velocity  diverges 
from  that  of  the  Gutenberg-Eullen  B  mcdel  and  tends  tcviard  the  phase 
velocity  curve  of  the  Bullen  A  model,  finally  meeting  the  Gun.’‘e  at 
about  550  seconds.  The  observed  curve  follows  the  theoretical  curve 
for  the  Gutenberg-Bullen  A  model,  at  first  slightly  below,  then  slightly 
above  ui^cil  350  seconds.  The  variations  in  this  period  range  frem  the 


theoretical  curve  are  orily  slight,  but  they  ,ire  consistent  aisd  do  not 
show  Sity  random  character.  Below  350  seconds  the  observed  phase  velocities 
diverge  doimwards  from  the  theoretical  phase  velocities.  Once  again  they 
approach  the  Gutenberg-Bullen  B  model  but  continue  right  through  tlas 
phase  velocity  curve  arvi  near  230  seconds  the  observed  ph^se  velocities 
are  lower  than  those  calculated  for  either  ir^del.  The  observed  phase 
velocities  in  the  period  rarige  of  200  to  liOO  seconds  are  in  excellent 
agi-eement  with  tte  values  obtained  from  surface  waves  by 

D.  Short  Period  Seismic  laves 

Thie  propagation  of  short-period  oceanic  surface  waves  vdth  pre¬ 
dominant  periods  between  5  arjd  20  seconds  was  studied  for  a  large  number 
of  pa'ihs  in  the  Atlantic,  Iridian  and  Pacific  oceans.  In  the  ocean  basins 
these  waves  are  controlled  largely  by  the  sedimentary  layer.  The  ampli¬ 
tudes  of  the  short  -period  waves  are  greatly  dependent  upon  the  riature  of 
the  propagation  path  ard  on  the  properties  of  the  source.  Surface  waves 
with  periods  between  5  and  10  seconds  are  often  observed  for  paths  that 
cross  the  contir^antal  margins  of  the  Pacific,  but  are  recorded  only 
rarely  for  similar  transmission  paths  in  the  Atlantic.  Probably  this 
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effect  is  d.^e  to  the  differenoes  in  the  stractural  coni  igi^rat ion  of  the 
margiriSj  particularly  to  such  differences  in  the  sedifscntary  strata* 

Seisniic  refraction  arid  reflection  noasurenients  indicate  that 
most  of  the  oceanic  paths  for  which  the  5  to  10  second  waves  are  observed 


are  characterised  by  sedimentary  thiclmesses  th^at  do  not  average  more 
than  a  few  tenths  of  a  kilometer*  The  prederdnant  periods  of  the  shert- 
neriod  wave  train  are  increased  to  as  much  as  12  to  20  seconds  for  paxhs 


that  traverse  the  thick  sediments  of  the  Argentine  Basin, 


inis  inc 


rease 


in  period,  which  vjbs  predicted  theoretically  for  regions  of  i*elatively 
thick  sediments  j  indicates  that  the  low-rigidity  sedi.ments  play  a  prominent 
role  in  determining  the  character  of  the  short-period  wave  train.  The 
increase  in  the  predominant  periods  cf  the  waves  associated  v?ith  the 
first  Love  and  first  shear  mode  also  accounts  for  the  absence  of  the  5 


to  10  second  waves  in  areas  of  thick  sediments  such  as  the  continental 


margiTiS  of  the  Atlantic, 

Dispersion  data  for  the  Rayleigh,  first  lo^re  and  first  shear 
modes  were  used  in  conjunction  with  reflection  and  refraction  results  to 
estimate  the  average  sh^ar  velocity  in  the  sediments  of  the  Argentine 
Basin,  The  average  shear  velocity  in  the  upper  0,5  km  of  sediments  is 
about  0,2  to  0,li  km/sec,  arid  the  velocity  in  the  kilonseter  of  sediments 
below  this  is  about  0,5  to  0,7  ion/sec,  SedimEntary  shear  velocities  of 
a  few  tenths  of  a  kilometer  per  second  were  also  obtained  for  paths 
alor,g  which  the  average  sedimentary  thicknesses  are  a  few  tenths  of  a 


kilomieter. 


■HHUlllll'’ 


For  group  velocities  between  ii«2  ar^i  3,h  km/sec,  the  particle 
motion  hi  the  first  shear  mode  is  retrograde;  ratios  of  horisontal-to- 
vertical  motion  as  large  as  3.5  are  observed  at  island  ard  coastal 
stations , 

Seismic  waves  with  periods  of  10  to  30  seconds  and  with  group 
velocities  of  h.O  to  h.li  km/sec  are  sometimes  recorded  on  tb^e  vertical 
components  of  long  period  seismographs.  These  v?ave3  may  be  explained  as 
higher  inodes  of  the  Rayleigh  type.  Oscillatory  waves  with  periods  of  5 
to  10  seconds  are  observed  to  follow  the  P  wave  in  many  oceanic  areas. 
These  arrivals  are  attributed  to  a  type  of  leaking  mode  that  results 
from  the  imiltiple  reflection  of  SV  waves  in  a  low- rigidity  layer, 

£♦  Epicenter  Relocation 

Seismicity  studies  were  carried  out  for  the  follcviing  areas* 
the  South  Pacific  Ocean;  the  Arctic  north  of  60°N;  the  G\ilf  of  Aden, 

Red  Sea  ard  Arabian  Sea;  the  African  Rift  Zone;  and  the  Central  Irdian 
Ocean, 

The  method  of  precise  epicentral  determinations  consists  of 
relocating  the  epicenters  of  earthquakes  using  a  computer  program,  A 
program  written  by  Bolt  (Geoptysical  Journal,  I960)  has  been  modified 
for  this  study,  rdre  than  200  P  and  PKP  readings  were  used  for  some  of 
the  larger  earthquakes  stidind.  All  known  earthquakes  within  the  above 
regions  were  relocated  for  the  period  1955  to  1962,  irith  the  exception 
of  the  South  Pacific  (195?  to  1963). 

Mien  previous  determinations  of  epicenters  were  plotted  on  a 
map  for  the  mid-oceanic  ridges,  it  was  fourd  that  the  data  scattered 
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over  a  b=elt  approxir^ately  200  to  500  km  In  width.  After  the  epicenters 
were  relocated,  the  data  were  concentrated  within  50  km  of  the  crests 
of  the  oceanic  ridges,  mch  information  concerning  the  geological 


structure j  the  position  of  the  mid-oceanic  ridge,  and  the  existence  of 
fracture  sones  can  be  learned  from  a  study  of  precise  locations  of  epi¬ 


centers.  In  each  of  the  areas  studied,  information  was  available  for 
more  eartliquakos  in  the  interval  1955  to  1$62  than  was  available  for 


the  previous  50  years. 


If  these  results  are  typical  of  seismic  belts  ard  of  small 


shocks  as  well  as  large  ones,  this  study  is  of  considerable  importance 


to  the  VELA  program. 


F»  Study  of  the  Deep  Argentine  Earthquake  of  December  8,  196^^ 
1.  Observed  travel  times  for  multiply-reflected  ScS  vraves. 
^faltiply-reflected  waves  from  the  December  8,  1962  deop-focus 
earthquake  in  Argentiria  vjere  recorded  very  clearly  at  stations  of  the 
Lament  IGY  and  USCGS  standardized  world -wide  networks,  ard  at  Gansdian 
stations.  The  data  gathered  from  this  earthquake  for  the  multiply- 
reflected  ScS  and  sScS  were  used  to  construct  ard  to  extend  to  shorter 


eoicentral  distances  the  travel  times  of  these  phases.  The  observations 
of  the  ScSScS,  ScSScSScS,  sScSScSScS  and  ScSScSScSScS  have  enabled  us  to 
obtain  the  travel  times  down  to  an  epicentral  distance  of  9^.  The 
existing  tables  so  far  shov;  the  travel  time  of  ScSScS  down  to  an  epi¬ 
central  distance  of  90°  only,  A  paper  on  this  subject  was'  rresehted 
at  the  Spring  196^  AGU  r^eting  in  Washington,  D.  C, 

In  the  process  of  plotting  the  ScSScS  travel  times,  a  mistate 
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was  found  in  the  ?i.iterature  for  the  ScSScS  travel  tinss  (Gutenberg 
and  Richter,  "Donnees  Relatives  a  L’ Etude  des  TrembleKients  de  Terre 
a  Foyer  ProfoM”,  Pablic^tions  du  Bureau  Central  Seismologique  Inter¬ 
national  1937,  and  in  Richter,  C»,  Elementary  Seismclo^,  Freeman  &  Co,, 
195B,  p,  68ii. 

2.  Rigidity  and  Q  factor  determination 

The  S,  ScS,  sScS,  ScSScS,  sScSScS,  ScSScSScS,  sScSScSScS  and 
ScSScSScSScS  have  been  recorded  very  clearly  at  some  of  the  IGi  and  uSCGS 
stations  arouTid  the  world,  and  at  Canadian  stations.  These  phases  have 
been  digitized  and,  using  a  Fourier  transform  program  for  the  IBM  1620, 
their  specura  have  been  obtained.  A  program  has  been  written  for  the 
same  computer  in  order  to  obtain  the  path  length  of  the  body  wave  under 
consideration.  This  program  uses  a  given  velocity  vs.  depth  distribution, 
and  coirpiutes  the  path  length.  A  program  which  calculates  the  reflection 
coefficients  at  the  core-mantle  boundary,  aixi  at  the  free-surface,  is 
used  in  order  to  correct  the  observed  spectrum  of  these  phases.  It  has 
been  found  in  the  process  of  computing  the  absorption  coefficient  that 
there  is  an  ejctra  term,  which  has  not  been  taken  into  account  by  Press, 
which  arises  from  the  geometry  of  the  problem.  For  an  ScS  the  geometrical 
factor  for  reflection  at  a  spherical  surface  is  divergent  in  nature,  and 
it  is  indepurdent  of  the  physical  properties  of  the  material  or  wave¬ 
length,  This  same  factor  of  reflection  at  a  spherical  surface  is  con¬ 
vergent  in  nature  for  the  ScSScS  and  higher  order  reflections  at  the 
core -mantle  boundary, 

fhiltiply-reflected  ScS  phases  (of  the  SV  type)  recorded  at 
close  stations  axe  used  in  the  calculation  of  the  rigidity  of  tte 


outer-core  cf  the  earth.  This  calculation  of  foe 
stations  at  larger  distances  is  under  stuct/* 

3.  Focal  mechanise  stvdy 


rl^iaitv  usir 


A  fault-plane  solution  lias  toen  mde  using  readings  of  coin- 
pression  and  dilatioris  of  36  stations  around  the  ‘^orld.  Dr.  Hodgson, 
of  the  Dominion  Oh=seru^atory,  has  run  a  fault-plane  solution  uith  a 
program  for  the  IBM  1620,  and  he  has  obtained  a  fault-plarse  solutxon 
very  similar  to  the  one  obtained  by  this  invastigator  for  the  same 
earthquake  of  December  8,  1962,  (a  single-couple  type  of  frechanism). 

P  atyi  S  waves  at  11  stations  u’cre  di^’-itized  and  Fourier  trans¬ 
forms  viere  obtained  for  these  phases.  These  spectra  of  F  and  S  waves 
have  been  plotted  and  will  he  coiipared  with,  spectra  obtained  for  P  and 
S  waves  from  theoretical  consideiations  applicable  to  a  double-couple 
mechanism.  The  theoretical  development  has  been  carried  out  by  Dr, 

T,  rIatuKoto,  of  Lament  Geological  Observatory,  Usin^  the  spectrum  of 
the  P  and  3  waves,  a  radiation  pattern  for  a  given  spectral  component  will 
be  plotted  versus  asi.mufn,  (Amplitude  (T)  vs.  Azimuth).  (Also,  an 
amplitude  (T)  versus  distance  plot  will  be  used  to  infer,  if  possible, 
how  sharp  the  core -mantle  discontinuity  is.) 


The  spectrum  of  the  P  and  3  waves  will  be  use'^,  also,  to 
ascertain  ariy  azimuthal  effect  on  the  spectral  structure  of  these  phases. 
Kasahmira  has  dons  a  study  along  this  line;  however,  he  did  not  have  ar^ 
conclusions  as  to  the  observed  difference  in  spectral  structure  of  the 
S  phases.  Eon  Menahemi  has  daiimed  that  this  difference  may  be  due  to  a 
moving  source. 


G,  The  Cliaracteristic  timbers  of  ths  Semi-Diurnal  Earth 
Tidal  Compcnent  for  various  Earth  ^fodels 

A  study  of  the  static  deformation  of  “"te  earth  due  to  tidal- 

generating  potentials  to  evaluate  the  possible  cause  for  the  discrepancy 


between  the  accmrailated  observational  data  of  earth  tides  for  the 

principal  semi-diurnal  tidal  constituent,  MU,  and  the  theoretical  values 

o 

has  been  made.  The  characteristic  numbers,  h,  k  and^  of  the  earth 
have  been  calculated  for  earth  models  consisting  of  v'elocities  due  to 
Jeffreys  or  Guttenberg  vrith  Bullen  A  or  Bullen  B  density  distribution. 


The  computations  were  made  using  a  modification  of  an  earlier  program 


for  an  IBM  7090  vrritten  for  calculating  the  free  oeriods  of  the  earth. 

It  is  found  that  the  characteristic  numbers  of  the  earth  are 
nearly  independent:  (1)  of  the  presence  of  the  so” id  inner  core  for  a 
wT.dc  range  of  the  possible  values  of  the  outer-core  rigidityj  (2)  of 
either  the  Jeffreys  or  the  Gutenberg  velocity  models j  i.e,  of  the 
presence  or  absence  of  a  low  velocity  channel  in  the  upper  mantle,  A  com¬ 
parison  of  the  theoretical  gravirrtetric  factors  Q  calculated  for  a  conti¬ 
nental  and  an  oceanic  earth  model  with  the  observed  gravimetric  factors 
of  the  5 end-diurnal  tidal  constituent  obtained  principally  during  the 


International  Geophysical  Year  indicates  that  the  differences  in  the 
observed  values  of  G  appear  to  be  primarily  due  to  the  effects  of  the 
ocean  and  only  secondarily'  to  differences  in  the  variation  of  the 
crustal  etructxires.  There  is  virtually  no  effect  due  to  the  preserrie 
of  the  solid  iriner  core, 

H,  PL  VJaves 

A  paper  entitled  ’’Dispersion  arai  Attenuation  of  Leaking  Modes" 
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by  S.  Su,  J.  Dorman 
describes  a  study  of 
higher  leaking  modes 


rd  Jack  Oliver  is  nearing  completiGn,  This  c=accr 
elastic  n-av's  corresponding  to  the  fundamental  arri 
in  the  continental  crust-ir^ntle  waveguide.  The 


approach  is  both  theoretical  and  experimental.  The  theory  is  based  on 
Haskell’s  (1962)  matrix  formulation  for  motion  at  the  surface  of  a 


layered  half-space  as  a  result  of  a  body  u’ave  incident  on  the  layers 
xromi  the  half -space.  A  method  is  developed  for  obtaining  dispersion 
curves  and  attenuation  for  waves  of  the  leaking  mode  type.  Theoretical 
curves  and  attenuation  computed  by  this  method  for  a  nwiiber  of  earth 


models  are  then  compared  with  dispersion  and  attentuation  data  of 
Oliver  (196li)  and  v?ith  some  nev?  data  for  certain  paths  in  South  America. 
Theoretical  values  cf  attenuation  for  simrle  earth  models  are  in  reasonable 
agreement  with  that  measured  for  PL  propagation  across  the  United  States 
by  Oliver  (156h).  Theoretical  curves  are  consistent  with  data  for  watres 


of  the  PL  type  and  of  the  shear-coupled  PL  typ=e.  In  particular,  it  is 
demonstrated  that  dispersion  data  for  the  sheai’-coupled  PL  v;aves  are  of 
special  value  for  a  more  refined  determination  of  earth  structure. 


ttirthermore,  certain  waves  which  appear  to  be  associated  with  at  least 
two  higher  leaking  modes  were  discovered.  The  identification  of  exteruied 


wave  trains  in  the  portion  of  the  seismograms  between  the  P  and  S  waves 


for  shocks  at  moderate  epicentral  distances  with  the  leaking  modes 
offers  the  hope  that  substantial  portions  cf  many  seismograms  can  be 
explained  in  this  rsanner. 


I.  Spectra  of  Body  V/aves 

The  nature  of  the  spectrum  of  a  body  ’wave  has  been  studied 
with  special  emphasis  on  p-eaks,  troughs  and  corresporwiing  periods. 
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Theoreticallyj  the  possible  causes  by  which  the  peaks  and 
trough.s  can  be  generated  in  a  spectrum  of  a  body  wave  are  as  follows: 


1.  The  Source  ffechanism  at  the  Origin, 


a.  According  to  the  distribution  of  an  initial  stress 
at  origin,  peak(s)  and  trough(s)  may  be  generated  in  the  spectrum  of  a 
radiated  body  wave.  For  the  excitation  caused  by  a  double»couple  t3rpie 
of  stress  distribution,  the  trough  period  in  the  spectrum  (Tt)  (the 


period  corresponding  to  the  minimum  amplitude)  is  approximate jy  given  by: 

Tt  ■=  rudius  or  spherical  origin/shear  u’cve  velocity  (l) 

fcr  the  P'wave.  The  spectra  generated  by  various  other  distributicns 

of  stress  are  also  being  investigated, 

b.  The  i^pothesis  of  a  moving  source  introduces  a  series 

of  peaks  and  troughs  in  a  spectrum  of  a  body  wave. 

The  periods  of  che  troughs  are  functions  of  the  azimuthal 
angle  Oo  aixi  the  rupture  =  velocity  are  given  by 


(2) 


where  b  is  the  extension  of  the  fault  and  N  is  an  integer. 

2.  Selective  Reflection  and  Refraction  at  Multilayered 
Discontinuities 

If  there  is  a  multilayered  medium,  selective  reflection  and 
refraction  occur  as  a  result  of  interference.  The  reflection  coefficient 
has  been  comouted  for  a  freo  surface# 


The  coiiiputed  reflection  coefficient  shews  a  series  of  minima 
in  the  amplitude.  The  corresponding  periods  are  functions  of  the 
incident  angle,  the  elastic  contents,  the  densities  of  the  layers  and  the 
thickriess  of  the  layers. 


?7 
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3.  Tho  Effect  of  a  TL>riB  Window 

For  a  given  tism  function  3’Ct),  the  effect  of  a  tiite  wirdow 


used  in  a  spectral  analysis  can  be  expressed  by 

-  2^  I  -/V  1  -O-iVa  (1)^  ' 

L>'  ' 


(3) 


where 


(  laS)  is  a  computed  spectrim 
(lO^  ^  Fcarier  transform  of  f'  (t) 

-  i-2  (  represents  the  effect  of  a  ti,me  v;indox7. 


If  the  time  windovr  ^  {t )  is  given  by 


Q  (t)  =  1 

0 


Ti‘5  t 


0  t<T3^ort>T2 


ih) 


—  L-  2  E-'^/is  expressed  by 

(i;  5  <^2-^ 


iUJ  jh  4  Ti) 


?: —  (5) 

Therefore  pseudo  peaks  arid  troughs  of  the  spectrum  may  be  introduced 
by  this  phenomena. 

ii.  The  Effect  of  Superposition  of  multiple  Arrivals 
If  the  seism.ogram  consists  of  two  (or  more)  arrivals  of  suc¬ 
ceeding  body  waves  with  a  time  separation  t.  (such  as  P  and  _P,  PKP, 

P  ^ 

and  PIiP2j  suvi  S  arid  sS),  these  waves  will  interfere  and  give  a  resultant 
airp)litude  modulated  vmve  train. 

The  spectrum  of  this  resultant  wave  train  'f'(t)  is  given  ty 


X  (uAj) 


)4(t)  +  a-f  (t-^); 


dt 


where  a  is  u  constant  (  1  a  1  ^  l). 


(6) 


For  a  1,  this  can  be  approximated 


o(uj)  i  2  cci  /2)  (7) 


Therefore 


(UJ)  has  minimum  amplitudes  at 


¥  (f  ’  (6) 

where  N  is  an  integer. 

The  initial  movements  of  the  Chilean  earthquake,  May  22,  I960, 
have  been  analyzed  and  compared  with  the  ti»oretical  aspects* 

No  evidence  of  a  moving  source  given  by  (2)  was  found,  but  at 
Rio  de  Jarteiro  and  at  Mt.  Tsukuba,  it  is  found  that  there  is  a  good 
agreeisent  between  observed  trough  periods  and  computed  trough  periods 
from  eq.  (8)  as  follows: 

l)  At  Rio  de  Janeiro  (A  "  30*0°) _ _ _ 


Observed  trough  periods  (sec„) 

120. 

37.3 

22.5 

17.3 

Computed  trough  period  from 
eq.  (8)  for  t  «  $?,$  sec. 

115. 

38.3 

23.0 

16.5 

N 

1 

2 

3 

h 

From  the  J.  B.  Tables,  the  time  interval  t  between  PP  and  P  is 


given  as 


t  «  $7.5  sec 


a)  At  m,  Tsukuba  {£k^  153.2^) 


29. 


Observed  trough  periois  (sec.) 

(20.8) 

15.9 

9.3 

7.2 

Computed  trough  period  from 
eq.  (8)  for  t  =  23*3  sec. 

ii6,6 

15.5 

9.3 

6.1 

N 

1 

0 

3 

h 

From  the  J.  B,  Tables,  the  time  interval  t  betvxeen  PKP2  and 


PKP^  is  given  as 


t  ®  21.9  sec. 


This  inplies  that  the  use  of  the  trough  periods  caused  inter¬ 
ference  gives  a  goodd]iie  on  hovj  to  separate  m>iltiple  arrivals  of  body 
waves  with  a  small  time  interval, 

J,  Surface  VJave  Inversion  Studies 

A  paper  entitled  "Seismic  Waves  and  Earth  Structure  in  the 
Canadian  Shield"  by  J,  Brune  and  J.  Dorman  was  published  in  the  Bulletin 
of  the  Seismological  Society  of  America.  For  this  paper  careful 
measurements  of  phase  velocities  in  the  Canadian  shield  have  been  made 
in  the  period  range  3  to  90  sec  for  Rayleigh  waves  and  12  to  60  sec 
for  Love  waves  by  phase  correlation  of  wave  trains.  The  continental 
Love  wave  phase  velocity  data  are  the  first  to  be  reported  in  the 
literature.  The  phase  and  group  velocities  are  higher  than  yet  found 
in  ar^  other  continental  area,  indicating  relatively  higher  shear 
velocities  in  the  crust  and  upper  mantle.  For  paths  in  the  Canadian 
shield  a  prominent  Lg  arrival  with  a  velocity  of  almost  3,65  km/sec  is 
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observed  am  an  Sn  arrival  is  recorded  clearly  to  distances  of  about 
iiOCC  iSTi  with  a  velocity  of  at-out  h,12  kia/sec. 

A  theoretical  layered  irndel  of  t}®  crust  aral  upper  inantle 
consistent  with  the  various  types  of  data  has  been  derived  by  an  in¬ 
version  method  employing  least-squares  curve -fitting  of  phase  velocity 
data.  This  model  has  a  three-layered  crust  35,2  km  thick  with  shear 
velocity  increasing  to  about  '^,85  km/sec  in  the  lower  crust.  The  upper 
mantle  has  a  nigh  speed  layer  with  shear  velocity  h,12  km/sec  dtarn  to 
about  115  km  belo\f  which  the  low  velocity  channel  has  a  shear  velocity 
of  about  ii,5  km/sec  down  to  a  depth  of  about  31e  km.  At  greater  depths 
the  shear  velocity  closely  follows  the  Gutenberg  model. 

Higher  mode  phase  velocity  dispersion  curves  comcuted  from 
the  theoretical  model  are  used  for  ccsuputing  theoretical  seismegrams 
for  the  Canadian  shield.  These  theoretical  seismograms  possess  maix,’  of 
the  features  of  the  observed  Ig  arrivals,  showing  that  Lg  can  be  explained 
by  nonnal  mode  wave  propagation  in  a  sifip^i,  layered  crust. 

This  paper  shows  that  new  nethcds  of  nsasuring  and  interpreting 
surface  wave  dispersion,  combiri!d  with  travel  time  data,  can  provide 
detaile  i  ard  reliable  information  on  shear  wave  velocity  distribution 
down  to  depths  of  a  few  hundred  kilometers  in  regions  of  horizontally 
uniform  structure, 

K.  Use  of  Partial  Derivatives  of  Surface  Wave  Phase  Velocities 
using  energy  equations,  which  are  equivalent  to  furdamental 
equations  and  bouidaiy  conditions  in  siirface  wave  problems,  one  can  obtain 
analytical  expr"*ssicns  for  partial  derivatives  of  phase  velocity  with 
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respect  to  physical  p=arair«ter  changes  within  the  earln.  These  derivatives 
can  be  calculated  by  knowing  only  the  eigenfunction  correspondir^  to  tfes 
phase  velocity.  fJumerical  examples  were  worked  out  for  proposed 

by  Brune  and  Dorman  for  the  Canadian  shield,  the  Jeffreys  model  and  the 
Gutenberg  model  for  crust-mantle  structure*  The  application  of  these 
results  to  numerical  inversion  of  surface  wave  data  has  been  demonstrated 
for  the  Canadian  shield. 

Item  It  -  Development  of  seisnac  instrumentation  vrlth  increased  sensitivity 
in  the  long  period  range  aid  with  improved  discrimination  against  noise 
in  the  dominant  microseism  band. 

A.  Seismo.neters  with  Disnlacement  Transducers 
One  iiathcd  of  increasLug  the  sensitivity  of  seismometers  in  the 
long  period  range  is  to  use  displacement  transducei^  instead  of  velocity 
transducers.  Therefore  the  lament  Geological  Observatory  set  itself  the 
task  of  developing  reliable  displacement  transducers  for  long  period 
seismographs.  This  has  been  accomplished. 

The  system  consists  of  the  following;  three -conoonent  long 
period  seismometers  v/ith  a  free  resonance  period  of  15  seconds  equipped 
viith  displacement  transducers  of  the  capacitance  type  and  associated 
amiplifiers  to  give  a  sensitivity  of  25  mdllivolts  per  micron  ground  dis- 
placemient  amplitude.  The  system;  provides  feedback  control  circuitry 
for  maintaining  center  aligrjnent  of  the  I.OP-g-period  seismometers. 

Gravity  aM  tilt  infonuaticn  are  provided  by  mnn Itcring  the  feedback 
signal  required  for  servo  centering. 
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The  basic  seismoisetcrs  used  are  ir^jdified  Press-Ewing  type 
iustruirents  manufactured  by  the  lehr^r  ar^d  Griffith  Division  of  IMited 
Elcctro-Dynanics .  The  associated  electronics  have  the  follci-;iiig  specifi¬ 
cations  ? 


Output  level;  +2.5  volts  DC  +  50  percent  with  a  zero 
to  five  (5)  volt  capability 
Output  sensitivity; 

Seismic  output;  .25  mv/mu 

2.5  mv/mu 
25  mv/mu 

depending  on  gain  setting 

Feedback  output;  ,232  volts/secund  of  arc  or  290 

volts /cm/s ec^ 

Long-term  stability  J 


Seismic  output: 
Tidal  output  i 


S%  of  full  scale 
,5^  of  full  scale 
LiTiCarity ;  10^  of  best  strai^t  line 

Squlvalent  input  noise:  25  jav  p-p  for  90^  of  the  time 

FreouencY  response t  (Filter  Passband ) 

Seismic  path;  l/60  cps  to  10  cps 

Feedback  path;  DC  to  1/2350  cps 

Attenuation  factors;  1,  10  and  1000  (nominal) 


Briefly,  the  principle  of  operation  is  as  follows:  The  cuter 
plates  of  the  balanced  capacitance  transducer  are  excited  ty  an 
oscillator  through  a  center- tapped  transformer,  in  push-pull  fashion. 
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Relative  motion  between  these  and  the  center  plate  of  the  transducer  is 
picked  up  as  an  amplitude -modulated,  suppressed-carrier  voltage,  amplified 
and  detected  in  the  Synchronous  Demodulator,  The  output  is  filtered  to 
supnress  high  frequency  components  ard  then  fed  through  an  attentuatcr 
to  a  high-pass  (T|.  -  60  sec)  filter  to  remove  the  DC  component.  It  is 
then  further  amplified  ty  the  Seismic  Amplifier. 

The  output  of  th«  Demodulator  is  also  applied  to  a  low-pass 
(Tj  =  2350  sec)  filter  whose  output  is  used  for  continuous  servo-centering 
of  the  transducer  through  a  coil-magnet  assembly.  This  assembly  is  also 
used  to  provide  the  proper  damping  for  the  seismometer .  Resistor 
determines  the  amount  of  this  positioning  feedback  and  resistor  R^j 
provides  the  damping.  The  feedback  signal  is  proportional  to  long-period 
motion  of  the  seismic  mass.  The  other  side  of  the  seismometer  coil  is 
connected  to  a  2,5  v  supply  source,  so  that  the  reference  level  is  estab¬ 
lished  at  mid-scale. 

In-phase  ard  l80-degree  cut-of -phase  components  of  the  feed¬ 
back  signal  are  fed  into  a  circuit  which  centers  the  vertical  instrument 
should  the  instrument  drift  outside  its  operating  range.  The  system  is 
calibrated  ly  applying  a  current  pulse  into  the  seismometer  coil  through 
an  attenuator  working  in  reverse  fashion  from  the  attentuatcr  in  the 
seismic  path. 

These  instruments,  xrlth  and  without  twin-T  band  rejection 
filters  in  the  seismic  output  have  teen  in  continuous  operation  at  Lament 
for  several  months.  Three -component  seismic  outputs  and  three -component 
tidal  outputs  are  recorded  routinely. 
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These  instmments  demonstrate  excellent  mechanical  and 
electrical  stability  and  minimize  the  need  for  the  electromagnetic 
transducers  stardard  on  many  seismometers  which  introduce  an  inherent 
6  db  octave  drop-off  in  the  response  to  long  period  signals.  The  use 
of  instrujtKnts  of  this  type  will  allow  long  period  seismic  signals  to 
be  recorded  more  easily  than  heretofore.  The  tidal  outputs  record 
solid  earth  tides  routinely  and  at  the  same  time  provide  an  extreme 
long  period  seismic  response  so  that  long  period  mantle  surface  waves 
may  also  be  recorded.  Although  these  systems  are  currently  more 
experisive  than  standard  seismoneters,  the  technique  will  prove  more  th.an 
worthjrhile, 

B,  Ifee  of  Surface  Wave  Rejection  Filters  to  Record  Mantle 
Waves  of  Low  Order 

Following  the  successful  use  of  galvanometer  band-rejection 
filters  to  reduce  microseismic  ra)ise,  the  saii^  technique  is  being  used 

to  filter  ordinary  surface  waves  arid  thus  obtain  clear  'Teccrdings  of 
early  mantle  wave  trains  from  large  earthquates.  The  long-period  waves 
are  usually  so  obscured  by  the  ordinary  surface  waves  that  it  is  diffi¬ 
cult  to  analyze  them.  For  several  recent  large  earthquakes,  clear 
records  of  mantle  waves  of  order  R1  (travelling  directly  from  earthiuaire 
to  station)  and  order  R2  (travelling  indirect  path)  have  been  recorded. 
The  records  are  ideal  data  for  the  stuc^  of  source  mechanisms  of  large 
earthquakes,  i.e,  the  determination  of  such  information  as  length  of 
fault  and  orientation  of  forces. 

At  present  the  system  consists  of  a  normal  30-100  long  {»riod 
Press-Ewing  seistiiOKraph  system  with  three  band-rejection  filter  galva- 
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noffiBters  in  series.  Each  filter  galvanoir^ter  is  preceded  by  a 
3000 -^'-shunting  resistor.  The  value  of  these  shunting  resistors  is 
critical  in  determining  the  shape  of  the  response  curve  for  the  system. 
For  larger  values  of  shunting  resistance,  the  system  is  more  tightly- 
coupled  ifith  greater  sensitivity  but  a  more  distorted  response  curve. 

The  periods  of  the  filter  galvanometers  are  set  at  15,  25 
arid  35  secords.  The  choice  of  periods  is  also  critical  in  determining 
the  response  of  the  system,  V/ith  the  present  system  the  filter  galva¬ 
nometers  occasionally  overload  (rebound  against  the  steps)  for  large 
earthquakes.  It  may  be  necessary  to  add  one  additional  filter  galva¬ 
nometer  to  correct  this  situation.  The  system  as  presently  operated 
has  been  calibrated  by  analysis  of  its  transient  response. 

Item  5  -  Conduct  studies  and  develop  techniques  in  the  use  of  specialised 
instrumentation  for  seismic  data  recording  and  analysis.  Included  would 
be  investigations  of  tape  recording  systems,  filtering,  digital  analysis 
arsi  analogue  computers  for  the  purpose  of  improving  detection  techniques 
and  methods  of  analysis  of  seismic  records, 

A,  Analogue  Techniques 

Analogue  computational  techniques  may  be  used  in  conjunction 
with  present  seismic  instrumentation  to  extend  their  capabilities.  Also, 
only  by  the  use  of  analogue  techniques  will  it  be  possible  to  process 
the  flood  of  data  which  will  arise  in  connection  with  the  needs  of 
Project  VELA.  Several  new  techniques  for  analysis  of  seismic  data  have 
been  developed.  These  include  the  following  items: 

1)  For  many  problems  involving  successive  filtering  to  look 


at  different  portions  of  the  frequency  spectrum,  p^ticular  slopes  on 
the  seisnsDgra^  response  cturves  are  required.  In  a«iiition,  inverse 
filters  are  required  to  shape  the  instrumental  response  or  to  remove 
the  instrumental  effect  entirely  from  the  seismogram.  The  building 
of  inverse  filter  elements  is  particularly  simple  using  analogue  coipputer 
ccupcnents  in  the  design  stages,  Ihe  use  of  these  inverse  filtering 
techniques  will  become  progressively  more  iiicortant  as  more  iragnetic  tape 
data  become  available, 

2)  Successive  filtering  of  long-period  seismic  signals  has 
been  shown  by  Sutton  and  Pomeroy  to  be  hi^ly  effective  in  removing 
microseisrlc  noise  and  in  the  selection  of  desired  portions  of  the  fre¬ 
quency  spectrum  for  subsequent  analysis.  This  technique  which  enhances 
the  signal  to  noise  ratio  allows,  in  seme  cases,  the  identification  of 
waves  otherwise  obscured  ly  higher  amplitude  waves  of  different  period. 

It  is  important  in  {^rforming  analyses  of  this  type  to  remove  the  filter 
phase  shift  characteristics  from  the  seismogram.  This  is  done  ty  playing 
the  seismic  signal  through  the  filter,  recording  the  filter  otitput  and 
playing  the  rerecorded  sigr.al  back  through  the  filter  backwards  in  time. 
This  results  in  a  flat  phase  spectrum  bit  the  amplitude  resporise  of  the 
filter  as  it  affects  the  signals  is  squared, 

3)  The  output  from  the  inverse  filtering  arrangement  described 
urder  1)  can  be  squared  and  integrated  witti  optional  filtering  to  provide 
an  output  proportional  to  u  ^dt  which  is  a  measure  of  the  kinetic 
energy  associated  with  the  ground  nation  at  a  particular  station.  The 
value  of  this  integral  can  be  sunaned  for  all  three  components  of  ground 


motion  to  give  a  irsasure  of  the  total  kinetic  Bmrgy  observed.  This 
number,  in  turn,  gives  some  measure  of  the  suarce  ener^’’  but  of  course 


does  not  take  into  account  the  variation  of  body  u^ave  energy  vri.th  angle 
of  incidence  (with  epicentral  distance)  or  variations  in  surface  wave 
energy  vdth  asixnuth.  The  study  of  the  ener©’'  accumulation  as  a  function 
of  tiiug  shows  vride  variation  on  the  three  coimoonents .  The  .measure  of  the 
total  emrey  accumulated  if  a  narrow  band  filter  is  used  gives  a  rough 
measure  of  the  nm^er  speetmm  of  the  time  series  on  the  magnetic  tape. 

U)  If  the  time  series  is  imultiplied  by  sine  and  cosine 
functions  ard  the  separate  results  are  inte gloated  to  form  the  values  of 
C(t)  Fourier  sine  and  cosine  coefficients  result, 

which,  if  squared,  added  and  the  square  roots  taken,  give  the  Fourier 
amplitude  coefficients.  These  coefficients  have  been  derived  for  several 
earthquakes  and  for  the  USSR  nuclear  shot  of  October  30,  I961.  The 
spectrum  0^  the  n\mJlear  explosion  as  recorded  at  Palisades  and  as  derived 
using  the  above  techniques,  is  remarkably  siiidlar  to  digital  spectra  of 
IS3R  explosions  of  the  1058  series  as  recorded  at  Palisades  ard  derived 


tising  digital  computers.  The  exception  is  that  there  is  evidence  of  some 
stronger  long-period  generation  lor  the  1^1  series, 

5)  If  the  azimuthal  angle  of  a  given  seismic  event  is  known, 
a  rotation  of  axes  analysis  may  be  performed  using  simple  analcue 
components  so  tbiut  fine  transverse  and  longitudinal  records  may  be 
obtained.  Once  this  analysis  has  been  perforjaed,  the  longitudinal  and 
vertical  comnonents  may  be  multiplied  together  to  separate  compressional 
type  particle  miotion  from  shear  type  particle  mi?tinn  in  most  cases. 
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Several  earthquakes  have  been  identified  using  this  technique  and 
several  unusual  features  have  becone  apparent  which  were  reported  by 
Rjirnroy  ai«i  Sutton  at  the  Seismological  Society  of  America  meetings  in 
Los  Angeles  in  April  1962. 

B.  Relative  Excitation  of  Long  and  Short  Period  Waves 

A  study  was  made  of  the  relative  excitation  of  long  period 
and  short  period  seismic  waves  by  earthquakes  and  explosions  in  the 
California-Nevada  region.  The  Richter  magnitude,  M,  ^:hich  is  based  on 
the  amplitudes  of  short  period  waves  recorded  at  sites  less  than  6C0  km 
from  the  epicenter  was  used  as  a  measure  of  short  period  waves.  A  new 
parameter,  AR  -  the  sum  of  the  areas  of  the  envelopes  of  the  surface 
waves  on  certain  three -component  long  period  instruments  -  was  defined 
and  used  as  a  measure  of  the  long  period  vraves.  The  surface  wave  para¬ 
meter  was  determiTied  from  long-period  seismograph  records  from  I^sadena, 
Berkeley,  Ruth,  Reno  and  Palisades.  The  four  West  Coast  stations  provided 
good  azimuthal  sampling  of  most  of  the  events  studied.  Surface  waves 
recorded  at  Palisades  from  the  Fallon-Stillwater,  Fairview  Ifeak-Dixie 
Valley  after  shock  sequernes  were  also  studied. 

There  is  a  wide  range  in  the  excitation  of  surface  waves  for 
a  given  Pdchter  magnitude,  both  for  explosions  and,  especially,  for 
earthquakes.  For  ary  given  earthquake  there  may  also  be  a  wide  range 
in  the  excitation  of  surface  waves  as  a  function  of  azii^tuth,  an  indication 
of  an  asymmetrical  radiation  pattern.  Generally  speaking,  the  relative 
excitation  of  surface  waves  by  explosions  is  much  smaller  than  by 
earthquakes  of  the  saaie  Richter  magnitude.  Most  of  the  earthquakes 
studied  generated  surface  waves  5  to  10  times  greater  than  the  maximum 
observed  for  explosions  of  the  same  Richter  magnitude. 
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A  ccmpariscn  of  the  relative  excitation  of  surface  waves  in 
different  regions  indicated  a  very  lar^  regional  variation.  Shocks 
off  the  coast  of  northern  California  generate  relatively  large  surface 
waves  whereas  shocks  near  Laguna  Salada  in  Baja,  California,  often 
generate  very  small  surface  waves. 

An  empirical  curve  was  derived  for  the  variation  of  the  surface 
wave  parameter  AH  as  a  function  of  distance.  It  was  sug^sted  th*at  the 
parameter  AR  may  he  used  to  estimate  the  surface  wave  magnitude  M  and  a 
formula  is  derived  for  doing  this. 

The  results  of  this  study  were  submitted  to  the  Journal  of 
Geophysical  Research  for  publication  under  the  title  ’’Relative  E^iita- 
tion  of  Surface  Waves  by  Earthquakes  and  liiderground  Explosions  in  the 
Galifornia-Nevada  Region”,  James  Brune,  Alvaro  Espinosa  and  Jack  Oliver, 
Prom  the  study  of  the  relative  excitation  of  long  period 
surface  waves  hy  earthquakes  and  explosions  in  the  Calif ornia-Nevada 
region  reported  above,  it  was  concluded  that,  for  a  given  Richter 
magnitude,  contained  xinderground  nuclear  explosions  generate  relatively 
small  surface  waves  compared  to  earthquakes.  It  is  suggested  that  it  may 
be  possible  to  use  this  fact  as  a  criterion  or  diagnostic  aid  to  identify 
certain  events  as  earthquakes.  It  roust,  of  course,  always  be  borne  in 
mind  that  the  results  described  here  arc  based  on  a  relatively  small 
sample  of  events  from  a  relatively  limited  number  of  regioriS.  There  is 
no  guarantee,  for  example,  that  a  nuclear  explosion  fired  under  different 
circumstances  will  not  generate  long  period  waves  even  though  such  a 
result  is  not  predicted  ty  elementary  theory. 
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As  an  example  of  hew  the  above  criterion  mi^it  be  used,  we 
have  iseule  a  statistical  st\xiy  of  the  number  of  earthquakes  studied 
which  produce  a  value  of  AR  1,  2,  5  and  10  times  the  maximum  expected 
for  undergrouiai  explosions  at  a  corresporaiing  distance  and  Richter 
magnitude.  For  a  given  natural  event  and  safety  factor,  e,g,  1,  2,  5 
or  10,  the  criterion  is  said  to  fail  (i.e,  fails  to  ideatiiy  the  event 
as  a  nat\3?al  event)  if  no  observed  value  of  AH  is  equal  to  or  greater 
than  the  maximum  value  of  AR  expected  for  explosions  at  that  distance 
and  magnitude,  multiplied  ty  the  appropriate  safety  factor.  Otherwise, 
it  is  said  to  succeed.  We  have  limited  the  magnitude  range  to 
3.91  M  5.5.  The  data  for  M  <  3.91  less  reliable  because 
of  the  spacing  of  the  stations  and  the  microseism  noise  level.  Of  the 
77  earthquakes  studied,  63  recorded  a  value  of  AR  greater  than  the  maxi¬ 
mum  expected  for  explosions. 

The  smaller  the  ratio  of  natural  events  not  excluded  to  the 
total  number  of  natural  events,  the  more  successful  is  the  method  in 
identifying  natural  events.  Depending  on  the  safety  factor  ard  magnitude 
range  chosen,  this  ratio  varies  from  iil/77  "  .53  for  a  factor  of  10  and 
magnitude  range  3.91  -  5.5  to  l/l6  »  ,06  for  a  factor  of  5  ard  magnitude 
range  U,71  -  5.5* 

If  the  lower  California  earthquakes  are  excluded  from  con¬ 
sideration,  in  the  magnitude  range  It. 31  -  5.5  the  ratio  is  0/21  ^  o. 
for  a  safety  factor  of  2,  1/21  «■  .C5  for  a  safety  factor  of  5,  and 
3/21  ®  .lit  for  a  safety  factor  of  10, 


The  reason  for  the  great  excitation  of  surface  v/aves  from 
the  earthouakes  off  northern  California,  and  the  very  low  excitation 
for  certain  earthquakes  in  lower  California  is  not  known,  but  could 
be  partly  due  to  different  source  nsechanisins ,  depth  effects,  poor 
magnitxide  determination  or  different  orientation  of  the  radiation 
pattern  with  respect  to  the  long  period  stations. 

Palisades  long  period  records  were  also  studied  for  earth¬ 
quakes  in  this  region.  Most  earthquakes  above  magnitude  htS  recorded 
long  period  surface  waves  and  fell  above  a  safety  factor  of  2,  and  many 
above  a  safety  factor  of  10. 

C,  Subsurface  Recordings  of  High-Frequency  Seismic  V/aves 

1,  Instrumentation 

In  late  1961  a  vertical  seismometer  with  electronic  amplification 
was  installed  at  the  Lairont  Seismograph  station  located  on  the  l85C-foot 
level  of  the  New  Jersey  Zinc  Company  mine  near  Ogdensburg,  New  Jersey, 

The  system  employs  a  Hall -Sears  K5-10  transducer  with  a  natural  fre¬ 
quency  of  two  cycles/sec,  an  SIE  GTR-200  amplifier  system,  a  high-speed 
photographic  drum  recorder,  and  a  two-channel  magnetic  tape  recorder. 

The  original  drum  recorder  operated  at  120  mm/minj  since  April  1^3> 
a  Sprengnether  auto-recorder  has  been  used  with  a  paper  speed  of  300  mm/mip. 
In  both  cases  the  SIE  amplifier  drovs  a  2CO  cycle/sec  "pencil"  galva¬ 
nometer,  Jfegmtic  tape  recording  of  the  SIE  amplifier  output  was  done 
on  a  two-channel  magnaohord  recorder,  with  Ifaemotron  pulse  frequency 
modulation  units  used  in  conjunction  to  pro^ride  a  pass  band  from  D,C, 
to  300  cycles/sec.  Recording  was  done  on  I/I4"  tape,  at  7  l/2  in/sec, 

A  ten- inch  reel  would  record  two  hours  of  data. 
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The  transducer-amplifier  'System  has  been  calibrated  In 
tteee  steps.  First,  a  irethod  employing  a  ’‘-illmore  bridge  yielded 
the  relative  response  of  the  transducer-aiaplifier  system.  The  methcxi 
takes  into  account  tl^  frequency  dejisndent  impedance  of  the  SIE  input 
transformers  at  frequencies  near  and  below  1  cycle/sec.  Second,  the 
electrodynamic  constant  of  the  transducer  was  measured  by  observing 
the  dajTping  as  a  function  of  the  resistance  terminating  tte  transducer 
ceil  winding.  Third,  the  gain  of  the  amplifier,  expressed  in  terms  of 
a  galvanometer  deflection  per  input  voltage,  was  measured  for  each 
record  by  means  of  a  test  signal  from  an  oscillator  within  the  SIE 
amplifier  system. 

The  over-all  magnification  of  the  system  with  photographic 
recording  increased  at  6  db/octave  in  the  frequency  barai  between  about 
3  to  100  cycles/sec .  The  system  was  U5t,ally  operated  at  a  magnification 
of  8  nilii-..-'!  at  20  cycles/sec.  During  the  quietest  recording  times,  the 
amplifier  noise  was  a  factor  of  ten  lower  than  the  equivalent  ground 
noise  over  the  freqiiencj  band  used, 

2,  Data  Analysis 

Ine  most  outstanding  characteristic  of  the  seismic  motion 
measured  at  the  mine  station  is  the  large  amplitude  disturbances  during 
daytime  and  early  evening  that  originate  from  milling  and  rock  crushing 
operations  associated  with  tiie  mine  and  with  a  limestone  quarry  nearby. 
The  disturbarces  f:om  these  sources  are  especially  predominant  at  fre¬ 
quencies  above  10  cps.  Because  of  -.his  high  daytime  noise  level,  irnst 
recording  was  doric  in  the  middle  of  the  night.  During  such  quiet  tines, 
approi'inately  800  hours  of  photographic  recording  and  about  20  hours 


of  tape  recording  were  obtained  on  the  l850-foct  le^l  with  the 
geophone -suuplifj^r  system.  In  addition,  about  20  hours  of  both  types 
of  recording  were  obtained  during  daylight  hours. 

Investigation  has  been,  so  far,  focused  on  the  quiet  night- 
tiiffi  periods.  Studies  of  the  continuous  seismic  backgrouTid  noise,  as 
well  as  the  various  discrete  signals  of  natural  and  c\dtural  origin, 
are  in  varying  states  of  progress.  Results  from  the  stidies  are 
presented  below, 

(a)  Continuous  seismic  bacliground  noise:  Several  two-minute 
samples  of  the  quiet  nighttiii©  bacl^ro’and  noise,  recorded  on  magnetic 
tape,  have  been  analyzed  with  one  octave  band  pass  filtering.  The 
resulting  amplitude  per  octave  spectra  show  a  sharp  drop-off  in 
amplitude  ty  a  factor  of  1C  from  0,5  cycles/sec  to  about  2  <ycles/sec. 
At  higher  frequencies  the  amplitxide  falls  off  at  a  rate  of  about  -6 
db/octave.  The  analysis  was  carried  ort  to  about  1^0  cycles/sec.  At 
higher  frequencies  the  noise  level  of  the  system  obscures  the  results. 

The  Lliiversity  of  Hichigan  recorded  and  analyzed  a  sample  of 
ground  noise  at  a  surface  site  very  nearly  above  the  1050-foot  level 
mine  site.  Their  results,  after  adjustment  to  differences  in  the 
analysis  procedur  and  collared  with  the  l850-foot  level  measurements, 
show  that  at  frequencies  less  than  oi®  cycle/sec  tte  amplitudes  on  the 
surface  and  at  depth  are  comparable,  but  that  at  hi^er  frequencies  the 
noise  level  unlergrcund  drops  off  faster,  reaching  a  level  one-fifth 
that  of  the  surface  noise  at  about  20  cycles/sec.  The  presence  or 
absence  of  the  two  to  three  cylces/sec  spectral  maximum,  observed  by 


ttii’VBrsity  of  Michigan  investigators  at  many  siarface  sites,  including 
Ogdensburg,  is  being  currently  examined  by  means  of  high  resolution 
analysis , 

(b)  Teles e is mie  body  phases;  The  only  body  phases  from 
distant  events  (  ^  )  v?hich  have  been  observed  to  contain  seismic 

energy  above  about  two  cycles/sec  originate  from  VJest  Indies  earthi^akes. 
Three  such  shocks  have  been  photographically  recorded  during  quiet  back- 
grouM  noise  tiires,  and  are  listed  below; 

1.  April  20,  1962.  N.  Coast  of  Haiti 
(USCGS  0  »  09:ii?:55.3,  20.6N,  72.2W) 

2.  Ifey  20,  1962,  N.  of  Puerto  Rico 

(15GGS  0  •=  I9t01;20,7,  20.9N,  66.0/) 

3.  March  13,  1963,  Dominican  Republic 
(ISCGS  0  =  10:39:19*1,  19.5N,  69.5W) 

The  peculiar  high  frequency  P  and  S  phases  often  observed  at 
eastern  U.S.A.  stations  from  such  shocks  were  well  recorded  in  all  of 
these  three  cases.  In  the  first  two  cases,  the  signals  were  not  vrell 
enough  resolved  on  the  photographj.c  paper  to  measure  the  frequencies 
in  the  signal.  Frequencies  as  high  as  10  cycles/sec  seemed  to  be 
present.  In  the  third  case  an  excellent  recording  was  obtained  at 
300  im/min.  Visual  inspection  showed  significant  ei^rgy  as  high  as 
5  cps  following  both  the  P  and  S  phases. 

Body  phases  from  other  distant  earthquakes  were  observed, 
but  with  a  frequency  content  no  different  than  that  iMicated 
recordings  on  more  conventional  short  period  seismotTaphn  such  as  the 
standard  station  Benioffs.  More  detailed  work  id.ll  be  done  in  order 


to  give  an  upper  ILndt  to  tte  propagated  erargy  with  frequencies  above 
3  cycles/sec  for  the  various  larger  earthquakes  that  occurred  Klien  the 
geophorie  was  recording. 

At  frequencies  arourKi  or^e  cycle /sec  thie  geophone-aiqplifier 
system  has  no  particular  advantage  over  conventional  cne-second  trans¬ 
ducer-galvanometer  systeiss.  During  daylight  hours,  as  discussed  above, 
the  noise  levels  at  higher  frequencies  are  high.  Consequently,  little 
has  been  done  with  regard  to  detecting  nuclear  shots  from  Nevada, 

(c)  Discrete  signals  of  local  or  regional  origin:  During 
the  quiet  nighttiine  periods,  several  types  of  signals  of  known  cultural 
origin  are  often  recorded.  Railroad  trains  numing  along  a  track  about 
ii  km  distaiit  at  the  nearest  point  produce  a  signal  lasting  several 
minutes,  which  is  marked  by  a  very  strong  three  to  four  cycles/sec 
component,  Tnis  dominant  frequency  is  also  observed  in  recordings  of 
trains  at  a  surface  site,  Sterling  Forest,  v^hich  is  about  7  km  from  the 
nearest  railroad  track.  This  phenomenon  is  currently  under  investigation; 
at  present,  it  is  believed  that  the  dominant  frequency  results  from  im¬ 
pulses  generated  bj?  the  train  wheels  traversing  the  track  joints. 

Quarry  blasts  are  occasionally  observed  during  the  earlier  parts  of  the 
night.  Host  quarry  blasting,  however,  is  dorie  during  daylight  hoars. 

Several  other  types  of  transient  signals  have  been  observed 
x^hcse  origins  at  present  are  not  definitely  known  but  are  suspected  to 
be  of  an  artificial  nature.  Very  short  duration  high-frequency  "blips” 
occxjT  occasionally,  about  once  a  night  on  the  average,  and  are  believed 
to  he  due  to  i^.^tall  local  readjustments  resulting  from  the  daily  mining 
activities.  Irregular,  rather  nondiscrete  disturbances  of  variable 


duration  and  frequency  content  also  occur  and  are  attributed  to  rock 
falls  within  the  nine  workinfs,  A  peculiar  pattern  of  disturbance, 
lasting  for  about  two  minutes,  arji  repeated  aLnost  exactly  several  tiines 
for  f!any  nights,  is  suspected  to  originate  from  a  rocket  testing  in¬ 
stallation  in  the  vicinity. 

Truck  traffic  on  a  major  highway  2  km  distant  produced  no 
perceptible  signals  on  the  l850-foot  level  instrument. 

The  records  h.ave  been  examined  closely  for  the  presence  of 
signals  that  ni^t  irjdicate  the  occurrence  of  small  earthquakes  in  the 
region.  733  hours  of  geophone  amplifier  recordings  during  tinass  of 
quiet  backgrourei  noise  corjiitions  were  searched  for  signals  with  the 
characteristic  P-S  signature.  Fourteen  events  were  fourvi.  On  the  basis 
of  investigations  to  date,  8  are  believed  to  be  small  earthquakes. 

Short  period  standard  station  records  obtained  in  the  mine,  as  well  as 
recordings  from  a  portable  Benioff  at  Sterling  Forest,  New  York  were 
also  examined  for  events  occurring  during  the  middle  of  the  ni^*tj  a 
total  of  2000  hours  for  the  two  instruments  yielded  6  events  suspected 
to  be  small  earthquairas ,  Pdchter  magnitudes  were  estimated  for  these 
l!i  events  and  locations  were  determined  for  those  large  enough  to  be 
recorded  at  several  stations.  Distances  to  the  events  from  Ogdensburg 
ranged  from  10  to  290  km,  and  magnitudes  ranged  from  -0.3  to  +2,9.  The 
located  epicenters  conform  with  historical  and  previous  instrumental 
determinations  of  areas  of  seismic  activity. 

Our  data,  combined  with  previousily  published  historical  and 
instrumental  data,  has  been  utilised  to  estimate  the  seismicity  of  the 
region  within  3CX)  km  of  0{|densburg,  New  Jersey,  If  the  data  is  plotted 


in  terins  of  a  number-  of  events  per  unit 


area  per  unit  tlnkj  with  magnitudes 
eQual  to  or  greater  than  seme  magnitudo  (the  iriiependent  variable), 
results  are  obtained  consistent  with  a  linear  relation  between  the 
logarithm  of  the  number  so  defined  and  the  .magnitude,  A  slope  of  -0.9, 
found  by  other  investigators  in  such  areas  of  high  seismic  activity  as 
Japan,  Southwest  Russia  aM  Southern  California,  seems  to  fit  the  data 
fairly  well. 

In  normalising  the  observed  number  of  events  to  a  number  per 
unit  area,  a  correction  has  been  applied,  peculiar  to  each  instrument, 
which  takes  into  account  the  variable  area  monitored  by  a  narticular 
instrument  for  a  given  magnitude  event.  The  correction  comes  from  an 
estimated  relation  between  the  miniimim  magnitude  event  detectable  arid  the 
epicentral  distance.  Recorded  events  Judged  near  the  threshold  of 
detecuability  were  used  to  estimate  the  relation,  for  the  most  part. 

The  geophone-amplifier  system,  having  the  greatest  magnification  in  the 
5-50  ops  frequency  band,  monitors  a  greater  area  for  a  given  magnitude 
event  than  the  conventional  Benioff  instruments  and  is  exceptionally 
sensitive  to  very  small  and  local  events.  Asada  in  Japan  and  Sanford 
and  Holmes  in  Hew  Ibxico,  both  using  instruments  similar  to  the  geo¬ 
phone-amplifier  systeni,  x^ere  able  to  detect  shocks  with  magnitudes  as 
low  as  "•2,  Their  results  were  employed  in  estimatirig  the  minimum 
detectable  magnitude  for  the  geo phone -amplifier  system  at  short  distances. 

The  resulting  estimate  of  the  frequency  distribution  of 
earthquake  occurrences  with  respect  to  magnitude  for  the  !^bw  York  region 
gave  6  X  10^  earthquakes  per  ye^  per  10^  km  with  magrd^  -ifs  equal  to  or 
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greater  than  zero.  Similar  estimates  of  the  seismicity  of  tvjo 
seisirdcally  active  regions  indicate  that  the  New  York  region  is  three 
orders  of  magnitude  less  active  than  the  fento  region  of  Japan  and  two 
orders  of  magnitude  less  active  than  Southern  California. 

Item  6  -  Conduct  tests  to  increase  the  knci/ledge  of  crustal  structure 
by  land  refraction  methods.  Field  studies,  using  a  portable  seismic 
recording  system,  will  be  conducted  at  selected  sites, 

Soir^  identifications  have  been  made  of  the  quarries  responsible 
for  blasts  recorded  in  connection  with  the  short  period  studies  in  a  deep 
mine,  previously  described  in  this  report,  for  the  purpose  of  future 
studies  along  these  lines.  Otherwise,  very  little  research  of  this  type 
has  been  carried  out  because  it  was  felt  that  research  in  other  directions 
was  more  suitable  fcr  the  aims  of  Project  VEIA, 

Item  7  -  At  a  site  to  be  designated  cy  the  Contracting  Officer,  establish 
a  long  period  seismograph  station  to  form,  tjith  the  existing  Palisades  and 
Wayneshurg  Stations,  a  tripartite  array  for  basic  studies  of  body  and  sur¬ 
face  waves.  An  additional  station  within  the  large  array  may  be  estab¬ 
lished  to  form  a  smaller  inner  array  if  these  studies  irdicate  it  to  be 
advantageous , 

No  suitable  site  was  found.  The  ^ditional  station  within 
the  proposed  array  was  set  up  at  Sterling  Forest  near  Tuxedo  Park,  Kev? 
York,  Since  this  item  was  first  proposed,  i/e,  as  well  as  others  working 
with  data  from  large  arrays,  have  experienced  diffi'^ulties  because 
horizontal  velocity  gradients  over  the  net  cause  variations  from  the 


expected  asiimiths  of  tte  arrivals. 


Accordingly,  it  is  likely  that  atr/ 


further  work  wi 
short  distances 


th  arrays  will  concentrate  on  arrays  composed  of  statioriS 
apart  and  located  as  much  as  possible  withdn  the  same 


geologic  province. 


Item  6  -  Develop  techniques  for  ar^alysing  microseismic  data  and  conduct 
studies  of  the  miicroseismic  spectra, 

A.  Direction  of  Approach  of  Micrcseisms 

Analogue  methods  have  been  used  to  determine  tl^  direction  of 
approach  of  Rayleigh  type  micrcseisms.  One-hour  samples  were  taken  for 
this  study.  The  azimuth  of  approach  determined  is  not  influenced  ap- 
pi'eciably  Love  waves  and  random  noise*  The  calculations  have  been 
perforinsd  on  a  TR-10  (EAI)  analogue  computer.  The  method  was  applied 
to  several  microseism  storms  recorded  at  Palisades,  The  results  are  in 
goed  agreement  with  those  of  earlier  studies,  but  they  are  far  more  pre¬ 
cise,  Differences  in  the  periods  of  microseisms  generated  by  different 
sources  allow  the  separation  of  multiple  source  problems  which  could  not 
be  handled  with  more  conventional  photographic  analyses .  Analysis  timiS 
is  significantly  reduced  by  playback  at  high  speed, 

B.  A  V/orldwidc  Storm  of  Microseisms  vdth  Periods  of  About 
27  Seconds 

On  June  6,  I96I,  between  about  1200  GOT  arid  2000  GOT,  a  storm 
of  microseismis  with  periods  of  about  25  to  2?  seconds  was  detected  by 
long  period  seismiOgraphs  at  marry  locations  tliroughout  the  world.  No 
storm  of  such  widespread  proportion  had  ever  before  been  reported.  Even 
though  the  sensitivity  and  the  geographical  distribution  of  long  period 
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instruicentation,  and  consequently  the  ability  to  detect  such  a  storiEj 
has  Improved  markedly  in  recent  years,  this  event  was  so  prominent  at 
some  stations  that  jt  seemed  unlikely  that  sirdlar  events  had  been  veiy 
frequent  in  the  reeeut  past.  There  may,  however,  have  been  storms  of 
similar  nature  but  smaller  size  which  have  gone  unnoticed,  for  many 
storms  of  long  period  microseisms  have  been  observed  at  individual 
stations. 

In  addition  to  its  widespread  distribution,  this  storm  was  un¬ 
usual  and  outstanding  in  certain  other  respects.  The  periods,  26  to  27 
seconds,  of  the  waves  were  considerably  longer  than  the  periods  of  micro¬ 
seisms  normally  observed,  even  on  the  few  seismographs  operating  with 
high  sensitivity  in  this  period  range.  On  such  instruments  the  back¬ 
ground  noise  usually  is  made  up  of  waves  with  periods  of  about  20  seconds 
or  less.  For  less  specialized  long  period  seismographs,  the  background 
noise  consists  of  ordiriaiy  microseisms  corresponding  to  the  well-known 
peak  of  the  earth  noise  spectrum  between  about  3  abd  10  seconds  (Brune 
and  Oliver,  1959). 

The  duration  of  the  storm  of  June  6,  only  a  few  hours  at  most 
stations,  was  much  less  than  the  duiation  of  most  microseismic  storms, 
Ultra-Sensitive  instruments  at  one  station  show  that  this  storm  actually 
had  a  considerably  longer  dxiration,  but  the  interval,  when  the  amplitudes 
were  very  large  was  very  brief.  During  this  interval  the  bandwidth  of 
the  frequency  spectrum  was  apparently  much  narrower  than  in  ordinary 
storms,  VJhereas  the  length  of  a  beat  in  an  ordinary  storm  is  usually 
5  to  10  cycles,  the  beats  of  the  June  6  storm  were  sometimes  as  long  as 
liO  or  50  cycles. 
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The  June  6  storm  also  had  the  property  that  the  periods  of 
the  waves  decreased  gradually  throughout  the  entire  duration.  At  most 
stations  vxhere  the  storm  was  observed  only  for  a  few  hours,  the  periods 
of  the  waves  decreased  from  about  27  seconds  to  about  26  seconds.  At 
Palisades,  where  ultra-sensitive  instruments  are  operated,  the  storm 
began  slightly  earlier  with  periods  of  about  28  seconds  arjd  continued 
through  the  interval  of  large  amplitudes  for  about  2  days  until  the 
periods  were  less  than  20  seconds.  The  gradual  decrease  of  period  with 
time  was  quite  certain  wherever  observations  were  good,  but  the  experi¬ 
mental  error  in  measuring  the  exact  period  at  aCiy  given  time  and  the 
slow  variation  of  peiiod  with  time  prevented  a  precise  determination  of 
the  velocity  of  propagation  cf  the  microseisms  between  distant  stations. 
However,  the  arrival  times  were  never  inconsistent  with  propagation 
bett^een  stations  as  seismic  waves. 

On  the  basis  of  particle  motion  studies,  the  seismic  waves 
appeared  to  be  predominantly  of  the  Rayleigh  type,  although  some  Love 
waves  might  have  been  present  as  well,  VJhen  the  waves  are  of  the  Ray¬ 
leigh  type,  the  direction  of  the  approach  could  be  determined,  and  ly 
combining  such  resxilts  from  several  stations  the  source  of  the  seismic 
waves  could  be  localized.  The  June  6  storm  appeared  to  have  originated 
in  or  near  the  southern  or  ec^atorial  Atlantic  Ocean,  possibly  in  the 
geraral  vicinity  of  the  Gulf  of  Guinea. 

Note  1,  In  the  performance  of  the  above  research,  the  con¬ 
tractor  will  continue  to  operate  the  widely  distributed  network  of 
long  and  intermediate  period  seismic  stations  that  were  initially  put 
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uiwo  operawiun  in  connection  vrlth  the  In-'^ernational  Geophysical  Year, 

with  such  changes  in  numbers,  location  or  instrumentation  as  may  be 

deemed  desirable  for  the  performance  of  the  contract. 

The  lon^i  period  seismic  stations  initially  put  into  operation 

during  the  IGY  were  continued  in  operation.  During  this  period,  the 

stations  at  Halle tt,  Antarctica}  Perth,  Australia}  Resolute,  North 

Canadian  Islands}  and  Suva,  Fiji  xfere  discontinued.  New  stations  were 

placed  in  operation  at  Reno,  Ifevada}  Toolangi,  Australia}  and  College, 

Alaska.  Two  sets  of  instruimsnts  will  be  placed  in  operation  at 

Wellington,  New  Zealand.  The  instruments  at  Agra,  India  x^ere  moved  to 
Delhi, 
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